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ABSTRACT 

We present computations for the accretion disk limit cycle model in an 
attempt to explain the empirical relation for dwarf nova outbursts between the 
peak visual absolute magnitude and orbital period found by Warner. For longer 
period systems one sees intrinsically brighter outbursts. This is accounted 
for in the limit cycle model by the scaling with radius of the critical surface 
density S max which triggers the dwarf nova outbursts. During the storage phase 
of the instability the accretion disk mass must be less than some maximum 
value, a value which scales with radius and therefore orbital period. When the 
instability is triggered and the accumulated mass is redistributed into a quasi 
steady state disk in outburst, the resultant peak optical flux from the disk is a 
measure of the total mass which was stored in quiescence. We compute light 
curves for a range in outer disk radius (or equivalently, orbital period), and 
find that our peak values of My are within < 1 mag of the observed relation 
M y (peak)= 5.64 - 0.259P orbita i(h) for 2 < P orbit ai(h)< 8. 

Subject Headings: accretion, accretion disks — instabilities — cataclysmic variables — 
stars: individual (SS Cygni) 
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1. INTRODUCTION 



Dwarf novae are a subclass of the cataclysmic variables — interacting binary stars 
with orbital periods of several hours in which a Roche lobe filling K or M dwarf secondary 
transfers matter at a rate M T into an accretion disk surrounding a white dwarf (WD) 
primary. Dwarf novae are characterized observationally by having outbursts of several 
magnitudes which recur on time scales of days to decades and which can last from ~ 1 day 
to several weeks (Warner 1995). The accretion disk limit cycle model has been developed 
and refined to account for the outbursts (Meyer & Meyer- Hofmeister 1981; for recent 
reviews see Cannizzo 1993a and Osaki 1996). In this model, material is stored up at large 
radii in a relatively inviscid accretion disk as neutral gas, and then dumped onto the central 
star as ionized gas when a certain critical surface density is achieved somewhere in the disk. 

Warner (1987) presented a thorough systematic study of the dwarf novae. He noted 
several interesting correlations between various attributes associated with the outbursts. Of 
particular interest is his finding of a relation between the peak absolute magnitude of dwarf 
nova outbursts and their orbital periods: My(peak)= 5.64 — 0.259P orb it a i(h)- In deriving 
this relation, Warner made corrections for distance and inclination for ~ 20 systems that 
were well enough studied to have reliable values for My (peak) and orbital period P rbitai- 
The distance determinations were made primarily on the basis of infrared fluxes of the 
secondary stars, as explained in Warner (1987). 

In this work we quantify the theoretical relation between My (peak) and P rbitai by 
running computations of our time dependent code which describes the accretion disk limit 
cycle model for a range in values of r outcr , the outer radius of the accretion disk. We 
compare our results to Warner's empirical relation. 
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2. BACKGROUND 

The physical basis for the limit cycle lies within the vertical structure of the accretion 
disk. In particular, the change in the functional form of the opacity at ~ 10 4 K coincident 
with the transition from neutral to ionized gas leads to a hysteresis relation between the 
effective temperature T e s and surface density £ in vertical structure computations carried 
out at a given annulus. This hysteresis leads to local maxima £ max and minima £ min 
in £ when the locus of solutions is plotted as T cff vs. £. The resultant series of steady 
state solutions resembles an " S" . Instability analyses of the equilibrium vertical structure 
solutions and their linearly perturbed states show that the upper and lower branches of the 
"S" curve, those portions where d log T eff /d log £ > 0, are viscously and thermally stable 
and can accommodate physically attainable states during quiescence and outburst. The 
"middle" part of the "S" where d log T e s/d log £ < is unstable and physically unattainable. 
During quiescence material accumulates at large radii in the disk because the viscosity is so 
low that there is little inward transport of the gas. The disk is far from steady state. When 
£ > £ m ax somewhere in the disk, a heating instability is triggered which begins an outburst; 
and when £ < £ m i n , a cooling transition front is started which controls the decay of the 
outburst from maximum light. As the disk goes from quiescence to outburst, the matter 
in the disk is redistributed from large to small radii, and the local M(r) profile becomes 
close to steady state, except for a strong outflow at large radii. The local maximum in £ 
determined from the vertical structure computations is 

£ max = 11.4 g cm' 2 rfo 05 mr ' 35 "Jf, (1) 

where r 10 is the radius in units of 10 10 cm, m x is the WD primary mass in solar units 
Mi/1Mq, and a co id is the value of the Shakura & Sunyaev (1973) viscosity parameter on 
the lower (or neutral) stable branch of the S— curve (Cannizzo & Wheeler 1984). (To avoid 
confusion, we use the subscript "max" to refer to properties of the accretion disk associated 
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with £ max , and "peak" to refer to either the peak flux or peak accretion rate during 
an outburst.) The values for £ max found by other workers are similar, reflecting slight 
differences in the equations of state, opacities, or treatments of the boundary conditions. 
The largest uncertainties may turn out to be due to the treatment of convection in the 
vertical structure. Within the confines of standard mixing length theory which is used in 
the vertical structure computations, convection may in fact be less important in dwarf 
novae accretion disks than was previously thought (Gu & Vishniac 1998). Nevertheless, 
even if convection is set to zero, one still obtains the hysteresis between T e g and S (see Fig. 
4 of Ludwig et al 1994). 

Cannizzo (1993b, hereafter C93b) commented briefly on Warner's finding of the 
empirical relation between the peak flux level in a dwarf nova outburst and the orbital 
period: My(peak)= 5.64 — 0.259P or bitai(h). C93b used the following argument to derive 
a scaling law for the rate of accretion onto the WD during the peak of a dwarf nova 
outburst. One can express the mass of the disk at the end of the quiescence interval during 
which time material accumulates in the disk as /M max , where M max = / 27rr<ir£ max is the 
"maximum mass" that the disk could have reached in quiescence if the disk were filled up 
to the level E max at every radius. Once the outburst has begun and progressed for a short 
time, the surface density profile adjusts from one which was non-steady state in quiescence 
to one which is in quasi-steady state in outburst. This basically involves a sloshing around 
of the gas from large to small radii, so that S qu i escence oc r and £ out b ur st oc r _1 . If this 
time interval is short, one may equate the mass of the disk at the end of the quiescence 
interval to the mass of the disk at the beginning of the outburst. Therefore if one integrates 
a Shakura-Sunyaev scaling for S ou tburst = S(r, a,M) over the disk as was done for the 
quiescent state by taking M out burst = / 27rr<ir£ ou tburst, sets this equal to /M max , and then 
inverts this expression to obtain M, one derives an approximation for the peak M value in 
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the disk during outburst 

Mpeak = 1.1 x 10~ 8 M e yr -1 (^p (j^f (Q , (2, 

(C93b), where ahot is the alpha value on the upper (or ionized) stable branch, and r oute r 
is the outer radius of the accretion disk (which is set by the orbital period). The values 
of the parameters entering into eqn. (2) have been scaled to the values which C93b 
found to be relevant for SS Cygni, a dwarf nova with P or bitai = 6.6 h. From Kepler's law 
-^orbital ° 3 > where a is the orbital separation. So if r outer oc a we basically have M pea k 
scaling as Porbitai) assuming the the other parameters in eqn. (2) do not vary with orbital 
period. If the visual flux were to scale linearly with M in the disk, then this would imply 
My (peak) oc —0.68 log P or bitab a different functionality than that observed. One additional 
consideration which comes into play for dwarf novae at increasingly longer orbital periods is 
that r oute r varies nonlinearly with a in the regime where the secondary star transitions from 
being less massive than the primary star to being of comparable mass. 

Warner (1995, see his Fig. 3.10) over-plotted eqn. (2) with the data taken from Warner 
(1987) to show that the analytical expression does a reasonable job of characterizing the 
observations. The analytical expression for M peak derived above is only useful in comparing 
with observations, however, if the variables appearing in eqn. (2) do not vary appreciably 
with radius. Of particular concern is the variable / which gives M outburst /M max at the time 
of burst onset. One might expect for / to vary significantly with orbital period -Porbitai or 
secondary mass transfer rate Mr, in which case the dependence on r outcr which appears in 
eqn. (2) would be misleading. Our aim in this work is to understand the function / by 
running time dependent models for a range in values of r outer and M T . 
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3. THE MODELS 



We use the computer model described in previous works (C93b, Cannizzo 1994, 
Cannizzo et al. 1995, hereafter CCL). This is a one dimensional time dependent numerical 
model which solves explicitly for the evolution of surface density and midplane temperature 
in the accretion disk. We carry this out by solving the mass and energy conservation 
equations written in cylindrical coordinates and averaged over disk thickness. The scalings 
which characterize the steady state relationship between T eS and £ were taken from the 
vertical structure calculations (Cannizzo & Wheeler 1984, C93b). For the WD mass we 
adopt Mi = 1M , for the inner disk radius we take r inncr = 5 x 10 8 cm, and for the number 
of grid points iV = 300. Our grid spacing is such that Ar = y/r. We compute the visual flux 
as described in C93b by assuming a face-on oriented disk and taking the standard Johnson 
V— band filters. C93b presents many tests of the numerical model to assess systematic 
effects. 

For the a parameter which characterizes the strength of viscous dissipation and 
angular momentum transport within the accretion disk we utilize the form given in CCL, 
a = a (h/r) n , where n ~ 1.5. CCL quantified the use of this form based on the observed 
ubiquitous exponential decays seen in soft X-ray transients and dwarf novae, but they 
were not the first to use it; it was introduced by Meyer & Meyer-Hofmeister (1984). 
The normalization constant a ~ 50 is based on the magnitude of the e— folding time 
constant, and the exponent n determines the functional form of the decay: n < 1.5 gives 
a faster-than-exponential decay, and n > 1.5 gives a slower-than-exponential decay. The 
reason for this particular functional form being the preferred one, based on a detailed 
examination of the departure from steady state conditions within the hot part of the 
accretion disk, has been recently provided by Vishniac & Wheeler (1996). There is 
one shortcoming in adopting this form for dwarf novae which must be addressed. This 
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shortcoming can be rectified if one imposes an upper limit on a, which seems physically 
reasonable. Figure 7 of CCL shows the failure of the form a oc (h/r) 1 - 5 to effect a change 
in the e— folding decay time in an outburst with outer disk radius, or equivalently, orbital 
period. For dwarf novae, however, it is well known that there exists a linear relation 
between the e— folding time constant associated with the decay of the dwarf nova outbursts 
and the orbital period (Bailey 1975). Therefore, the CCL a value cannot be applicable to 
dwarf novae without modification. Recent computations applying the form a = ao(hfr) 1 - 5 
to integrations of the vertical structure show that, if a is computed self-consistently within 
this formalism for dwarf novae, one derives large values of a in the "outburst state" — 
meaning the point on the upper branch of the log T eff — log £ curve which lies at the 
vertical extrapolation of the S max value (Gu & Vishniac 1998). In fact, the values are 
considerably larger than can be tolerated within the framework derived from the theoretical 
constraint imposed by the Bailey relation (Smak 1984, Cannizzo 1994) — this constraint 
being that «hot cannot exceed ~ 0.1 — 0.2. This limit will be reached, unfortunately, when 
h/r ~ 0.016 — 0.025 if we strictly take a = a (h/r) l - b . Since this value of h/r is exceeded 
on the upper stable branch of steady state solutions for dwarf novae, we conclude that the 
physical mechanism responsible for generating the viscous dissipation and transporting 
angular momentum in accretion disks must saturate to some au m it ~ 0.1 — 0.2 so as to give 
the Bailey relation (Smak 1984). Disks tend to be thinner in systems with larger central 
masses, generally speaking, so systems such as X-ray novae (and AGN) do not run into this 
limit (Gu & Vishniac 1998). This explains why X-ray novae do not show a "Bailey relation" 
(see Fig. 3 of Tanaka & Shibazaki 1996). In light of these considerations, we utilize the 
CCL a form in our computations, but in the limit of large a we do not allow it to exceed 
0.2. 

The fact that we see a system as a dwarf nova means that the mass transfer rate into 
the outer accretion disk from the secondary star cannot be too great or else the disk would 
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be in permanent outburst (Smak 1983). This reasoning must apply to the systems used 
by Warner (1987) in his compilation. Shafter (1992) considered the relative frequency of 
dwarf novae as a fraction of all cataclysmic variables in different period bins longward 
of the 2 — 3 h period gap in an attempt to understand the variation of the rate of mass 
transfer M T from the secondary star (feeding into the outer accretion disk) with orbital 
period which we observe in the dwarf novae. Shafter concluded that, by restricting our 
attention solely to dwarf novae, we are probably not sampling the long term Mr(-Porbitai) 
value which characterizes the secular evolution of cataclysmic variables as a whole (Kolb 
1993). Therefore we need not feel reticent in adopting an Mr(-Porbitai) law which serves only 
to ensure M T < M crit , where M crit is the value of the secondary mass transfer M T which 
must be exceeded for the entire disk to be stable in the high state of the "S" curve. Stated 
another way, our concern in this work is to adopt a law for Mt (-Porbitai) which will produce 
dwarf nova outbursts. This law is not related to that which characterizes the entire class 
of cataclysmic variables. From a practical standpoint, we find that our computed values of 
My (peak) are relatively insensitive to the specific M T at a given orbital period. 

We run computations for five values of r outer /10 10 cm — 2, 3, 4, 5, and 6 — while at the 
same time scaling the value of the rate of mass transfer M T feeding into the outer disk so 
that M T < M crit = M(E min ) ~ 10 16 g s" 1 (r outcr /10 10 cm) M m^ 87 (Cannizzo & Wheeler 
1984). For our canonical SS Cyg model, we take r outer = 4 x 10 10 cm and M T = 6.3 x 10 16 
g s _1 (C93b, Cannizzo 1996). At this value, the system is about a factor of 4 below M crit . 
Therefore, to determine M T for other r outer values we scale the SS Cyg M T by r^ ter . We 
also ran a second series of computations with the normalization on a factor of two 
smaller. 

Figure 1 shows the computations from the first series of models. In Figure la we give 
the light curves, and in Figure lb we give the accretion disk masses, shown in units of M max 
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for the relevant r outcr value. The assumption implicit in eqn. (2) that / is constant with 
orbital period appears to be quite good: there is no noticeable drifting of the disk mass 
relative to M max as r outer changes. Figure 2 shows the computations from the second series. 
The results are similar to the first series. The equilibrium disk masses shift to slightly lower 
values, as do the My (peak) values. 

Figure 3 shows the results taken from the previous figures of My (peak) versus orbital 
period. The results from the first series are indicated by the hatched area with hatched 
lines inclined at ±45° with respect to the whereas those from the second series 

have hatched lines inclined at 90° and 180°. The conversion between r outer and P or bitai was 
carried out using the fitting formula given in Eggleton (1983), assuming a secondary mass 
M 2 = O.lM P orbital (h). The hatched area accompanying each series shows the range allowed 
by taking the disk to fill between 0.7 and 0.8 of the Roche lobe of the primary. 



4. DISCUSSION AND CONCLUSION 



We have run models using the accretion disk limit cycle model for dwarf 
novae in an attempt to understand Warner's relation for dwarf nova outbursts 
My(peak)= 5.64 — 0.259P or bitai(h). As noted in C93b, the observed upper limit is a natural 
consequence of the "maximum mass" of the accretion disk that is allowed by the critical 
surface density S max . This value increases steeply with orbital period, therefore the amount 
of fuel available in a dwarf nova outburst also scales with orbital period. An important 
finding is that, for a constant value of M T /M crit , the value / which is the ratio of the 
accretion disk mass at outburst onset to the "maximum mass" is relatively constant with 
orbital period. This gives some confidence in the analytical estimate given by C93b. The 
fact that the theoretical variation of My (peak) with orbital period is flatter than might 
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have been expected from C93b's scaling is due in part to the variation of the orbital period 
with outer disk radius in the limit where the mass of the secondary star starts to become 
comparable to the primary mass. This effect contributes to the arc-like shape of the shaded 
regions shown in Fig. 3. 

The mean level of My (peak) in our models exceeds Warner's empirical line by 
~ 0.3 — 0.5 mag over most of the range shown. It is probable that our method for computing 
the V band flux is too crude to expect consistency with observations at this level of detail — 
for instance we do not include limb darkening in the models, an effect mentioned by Warner 
(1987). Also, we have utilized Planckian flux distributions for the disk. Wade (1988) has 
shown that, to varying degrees, both Planckian distributions and stellar or Kurucz type 
distributions fail to represent adequately the observations. Although Wade found that the 
Planckian distributions can account for both UV flux and UV color in a sample of nova-like 
variables, the failure of the model in other respects leads to the conclusion that "one 
cannot rely on model fitting to give the correct luminosities or mass-transfer rates (within 
an order of magnitude)." Unless there are systematic effects which depend strongly on 
orbital period, however, our computed My (peak)— P orbita i slope should have some physical 
relevance. When one goes beyond this to considering the normalization level of My (peak), 
it would seem that a better flux model must be utilized. 

It is also interesting to note the relative insensitivity of My (peak) to M T at a given 
orbital period. This would imply that the observed scatter at a given orbital period must be 
due largely to variations in distance and inclination. We varied the normalization constant 
on Mr by a factor of two between our two series in this work, and found the resulting 
My(peak) values to differ by < 0.3 mag. 

The soft X-ray transients — interacting binary stars containing a neutron star or 
black hole as the accreting object — also seem to obey a relation between My (peak) and 
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-Porbitai for outbursts. Van Paradijs & McClintock (1994) noted a relation between My and 
(Lx/^Edd) 1 ^ 2 -Porbitai' where L x /L Edd is the X-ray luminosity in units of the Eddington 
luminosity. Irradiation is a complicating factor in these systems, and by comparing the 
My values between Figure 1 of Warner (1987) and Figure 2 of van Paradijs & McClintock 
(1994), one can see that the X-ray binaries are ~ 1 — 7 mag brighter than the dwarf novae 
at maximum light. Clearly most of the optical flux is reprocessed X-ray radiation coming 
from large radii in the disk. For the black hole X-ray binaries, however, a large fraction 
of the difference comes from having larger disks due to the larger orbital separations, at a 
given orbital period (Cannizzo 1998). For the neutron star systems (in which the primary 
mass is ~ 1M© as in dwarf novae), the difference is entirely due to the irradiation. 
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and Ben Zellner at the Space Telescope Science Institute. JKC was supported through 
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FIGURE CAPTIONS 

Figure 1. (a) The light curves from the first series for which (top to bottom) the 
values (r outer /10 10 cm, M T /10 16 g s" 1 ) = (2, 1.04), (3, 2.98), (4, 6.30), (5, 11.3), and (6, 
18.1). The dashed lines indicate the approximate absolute magnitude of the secondary star 
My = 22.0 — 17.46 log P or bit a i(h) (Patterson 1984). The two values shown in each panel are 
consistent with a disk which fills 0.7 and 0.8, respectively, of the Roche lobe of the primary. 
In the first two panels these lines are off-scale toward fainter magnitudes; in the third panel 
only the upper line is visible. 

(b) The accretion disk masses corresponding to the light curves in (a), scaled to the 
"maximum mass" values (top to bottom) M max = 2.713 x 10 23 g, 8.013 x 10 23 g, 1.727 x 10 24 
g, 3.135 x 10 24 g, and 5.103 x 10 24 g. 

Figure 2. (a) The light curves from the second series for which the Mt values are a 
factor of two smaller than in Fig. 1. The r oute r values are the same. 

(b) The accretion disk masses corresponding to the light curves in (a), scaled to the 
same maximum mass values as in Fig. 1. 

Figure 3. The peak absolute magnitudes of the outbursts shown in Figs, la and 2a. 
The upper hatched region is for the first series and the lower hatched region is for the 
second series. The normalization on the M T values used for the two series is different 
by a factor of two, and yet the offset between the two regions is quite small. To convert 
between r oute r and Porbitai, we use the fitting formula for the Roche radius given by Eggleton 
(1983), assuming either 0.8 (left boundary of each hatched region) or 0.7 (right boundary 
of each hatched region) for the fraction of the Roche lobe which r outcr occupies. We also 
take M 2 = O.lM P orbital (h). The dashed line indicates My(peak)= 5.64 — 0.259P or bitai(h) 
(Warner 1987). Warner quotes an rms scatter about this relation of ±0.23 mag. 



